Twists and turns for metamaterials 
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We propose and verify experimentally a new concept for achieving strong nonlinear coupling be- 
tween the electromagnetic and elastic properties in metamaterials. This coupling is provided through 
a novel degree of freedom in metamaterial design: internal rotation within structural elements. Our 
meta-atoms have high sensitivity to electromagnetic wave power, and the elastic and electromag- 
netic properties can be independently designed to optimise the response. We demonstrate a rich 
range of nonlinear phenomena including self-tuning and bistability, and provide a comprehensive 
experimental demonstration of the predicted effects. 



Metamaterials research has grown rapidly over the past 
decade, exhibiting a wide variety of new wave phenomena 
U [H . Being, initially conceived in the domain of electro- 
magnetics the metamaterial concept also proved 
to be fruitful in other areas of physics [6|-l8|. Until re- 
cently, however, direct interplay between different types 
of physical effects within the same metamaterial was not 
considered, although mechanical control over electromag- 
netic metamaterial properties was employed in structural 
tuning [Si, llQ|. 

It turns out that introducing a mechanical degree of 
freedom into electromagnetic metamaterials leads to an 
interesting range of nonlinear effects, giving rise to a new 
class of magnetoelastic metamaterials [11] and to wide- 
band operation [l2| . The range of possible effects achiev- 
able in this way promises to be richer than in the promi- 
nent area of optomechanics [Hj , because the greater flex- 
ibility in metamaterial design overcomes the limits of 
available material functionalities, and offers wider pos- 
sibilities for optimisation. At the same time, the imple- 
mentation of magnetoelastic metamaterials PJJ remains 
challenging and in some cases, such as the conformational 
nonlinearity in resonant spirals 14], remains inaccessible 
for optics. The reason for this is that the magnetic forces, 
employed in the initial designs, are relatively weak, so 
such materials require either high power or extremely 
small elastic restoring forces, which poses a considerable 
manufacturing challenge. 

We recall, however, that earlier research on struc- 
turally tunable metamaterials ^ indicated that near- 
field interaction may significantly improve the tunability 
range, leading to various effects associated with near-field 
coupling [15j. In particular, changing the mutual orien- 
tation between the neighbouring elements has a profound 
effect on the resonator coupling and the structure of their 
modes [ll,[l7|. 

We therefore expect that the most efficient approach 
to implement dynamic coupling between electromagnetic 
and mechanical effects should rely on near-field interac- 
tion, which can have a powerful influence even for subtle 



changes in the mutual orientation. Instead of a consid- 
erable displacement of an entire array [9], it is sufficient 
to move the crucial parts of the resonant particles with 
respect to each other — for example, the gaps of the two 
coupled split-ring resonators — which involves a more 
gentle geometric alteration. 

Let us now introduce our novel concept: nonlinear 
metamaterials with intrinsic rotation. As a building ele- 
ment of the structure (see Fig. [I]), we consider two coax- 
ial split ring resonators (SRRs) with elastic feedback be- 
tween them. The rings are allowed to rotate about the 
common axis, while the elastic feedback is provided by 
connecting a thin elastic wire. Indeed, the use of elastic 
wires has a prominent history in physics, being utilised 
in the milestone achievement of the experimental demon- 
stration of light pressure by P. N. Lebedev [18]. Here, we 
employ the electromagnetic (EM) torque to construct a 
"light- driven" meta-atom, which enables us to modulate 
the resonant frequencies by twisting the rotatable ele- 
ment directly with EM waves. The component of the EM 
forces which twists the rings with respect to each other 





FIG. 1: Conceptual layout of a new metamaterial and its 
rotational "meta-atom". Incident wave propagates along y 
direction, having a linear polarisation with the electric field 
along x and magnetic field along z. The induced electromag- 
netic torque between the resonators changes the mutual twist 
angle between the rings, connected by an elastic wire. 
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is normally not the strongest among the forces involved, 
but the prominent advantage of using EM torque instead 
of collinear EM force to drive a meta-atom is that the ef- 
fective lever arm of the azimuthal EM force can be much 
larger than that of the azimuthal restoring force from a 
thin wire; this can effectively magnify the deformation in 
azimuthal direction by orders of magnitude. 

Technically, there are a number of ways to implement 
this general scheme; in our design, one of the rings is 
fixed to a substrate and the other one is suspended on 
the long wire. The three symmetrically positioned wires 
which attach the suspended ring to the string provide 
stability against tilt. In this design, therefore, the only 
favourable movement is the rotation of the suspended 
ring with respect to the fixed ring over the common axis, 
and all other mechanical degrees of freedom can be ne- 
glected. 

Suppose the initial position is such that the ring slits 
have a certain angle between them with respect to the 
common axis (Fig. [1]). An EM wave then induces a cer- 
tain distribution of charges and currents in the two res- 
onators, and the resonance of the system is determined 
by their mutual orientation [l6|. These charges and cur- 
rents also result in EM torque between the two rings [HI , 
which drives the suspended ring to rotate until the EM 
torque is compensated by the elasticity of the twisted 
wire. Meanwhile, the entire pattern of charges and cur- 
rents gets modified and the torque also changes, so the 
final stable equilibrium is only achieved via a complex 
nonlinear feedback. The twisted wire provides a restoring 
torque to balance the EM torque, so that we can control 
the twist angle (and thus the resonance) by changing the 
external field signal. 

An additional feature of the proposed design is that 
the nonlinear dynamics of the rotatable particle not only 
depend on the parameters of the wire, but also on the 
EM mode initially excited in the resonators, which is de- 
termined by the starting angle between the gaps. The 
latter can be made arbitrary, and an implementation has 
the possibility to deliberately adjust it, introducing tun- 
ability to the system. 

The above design, therefore, offers a tunable resonant 
nonlinear system with elastic feedback and, as we show 
below, yields a rich pattern of nonlinear response includ- 
ing self-tuning and nonlinear bistability, which are much 
stronger than those provided by using nonlinear semicon- 
ductor components. Below, we present a detailed theo- 
retical analysis of rotational meta- atoms. We then pro- 
ceed to the experimental results obtained with a fabri- 
cated prototype of the rotational "meta-atom" placed in 
a rectangular waveguide, and confirm all important fea- 
tures predicted by theory. Finally, we perform full-wave 
numerical simulations of the array of "meta-atoms" (i.e., 
metamaterial), and demonstrate that all the nonlinear ef- 
fects observed for the single meta-atom in the waveguide 
are qualitatively the same in the array. 



To start with, we use a semi-analytical model to study 
the dynamics of an isolated meta-atom in free space. It 
will subsequently be demonstrated that this model ex- 
plains all qualitative features of an array. As shown in 
Fig. [TJ the two coaxial identical SRRs are offset by a dis- 
tance s in the z direction, and the twist angle between 
them is 0. The incident wave propagates along the y 
direction, with its magnetic field in the z direction and 
electric field in the x direction. To study the nonlinear be- 
haviour of the rotatable meta-atom, we utilise an efficient 
analytical model based on the single mode approximation 
and the near-field interaction [16]. This model can pro- 
vide a reasonable prediction of the EM response as well as 
the optomechanical properties of the structures [HI . As 
our previous studies showed [loj, the current and charge 
of the SRR can be separated into frequency-dependent 
mode amplitudes and spatially-dependent distributions: 
J(r,o;) = —jujQ(uj)j(r) and p(r,(j) — Q((j)q(r). The 
mode amplitudes Qi^ can be obtained after solving the 
coupled equations: 

Qi = (£2F m - £iF s )/(F s 2 - Fl) 
Qi = (£iF m - £ 2 F S )/(F* - Fl), 

where E\ and £2 correspond to the effective voltage ap- 
plied to the lower (index 1) and top (index 2) SRRs by 
the external fields; F s and F m are the self and mutual 
impedance terms, with their explicit forms given in the 
supplemental materials [21] . 

Once the frequency-dependent mode amplitudes are 
known, we can calculate the EM torque experienced by 
the SRRs. Here, we are particularly interested in the 
torque on the top rotating ring: Mem = f v p( r 2) 1*2 x 
E+r2 x [J(r2) xB]d^2, where the integration is performed 
over the volume V2 of the top SRR. We decompose the 
total torque into two parts: external t orq ue M ext con- 
tributed by the external incident fields [19| , and internal 
torque Mi nt due to the near-field interaction between the 
two SRRs. We assume that the lower SRR is fixed while 
the top SRR is only allowed to rotate about the z axis. 
See the supplemental materials [2l[ for explicit expres- 
sions for the internal and external torque. 

Figures [21(a) and EJb) depict the mode amplitude Q2 
and the total EM torque Mem = M ext + Mi nt experi- 
enced by the top SRR as functions of frequency and twist 
angle 0. Since radiation losses are taken into account in 
the model [20] , we are able to accurately describe the evo- 
lution of the mode amplitudes, phases and lineshapes of 
the resonances. As expected, this chiral meta-atom sup- 
ports two hybrid resonances, which can be characterised 
as symmetric (lower frequency branch) and antisymmet- 
ric (higher frequency branch) modes, according to the 
symmetry of the H z component [16]. 

The directions of the EM torque at these two res- 
onances are also opposite. For the symmetric mode, 
= 0° corresponds to the configuration of highest po- 
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tential energy (unstable point), and thus the two repel 
each other once > 0°, until they come to the stable 
state at = 180°, while the reverse is true for the an- 
tisymmetric mode [22|. The evaluated external torque 
is about one order of magnitude smaller than the in- 
ternal torque, and the total torque is of the order of 
10 _10 Nm when the structure is pumped with a power 
density Pi=lmW/mm 2 . This is confirmed by the full- 
wave simulation (CST Microwave Studio) followed by 
calculation based on the Maxwell stress tensor, which 
yields the EM torque through a surface integral of the 
field components around the object [23| . 

The overall mechanism of achieving a nonlinear ef- 
fect is presented in Figs. [2fc,d). As an example, we 
choose a pump frequency (3.5 GHz) at the symmetric 
mode [regime denoted by the black dashed line in Fig. [2] 
(b)]. It can be seen that Mem is a Lorentz-like func- 
tion of the twist angle, while the restoring torques Mr 
under different initial twist angles #0 are approximated 
by linear functions (Hooke's law). The intersections of 
these two functions, Mem^^i) + Mr(^) = 0, corre- 
spond to the equilibrium angles e . However, only the 
angles with ^ [M EM (# e ) + M R (0 e )] < are stable. As 
the pump power P\ increases from zero to maximum and 
then reduces, the stable angles also change accordingly. 
With this method, we can numerically find a sequence of 
stable angles under different pump power P\. 

Since the EM torque is a nonlinear function of twist an- 
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FIG. 2: The principle of nonlinear response in rotatable 
meta-atoms. (a) The mode amplitude Q2 and (b) the EM 
torque Mem of the top rotatable ring, (c) The EM torque 
at 3.5 GHz for different pump powers from to 1 mW/mm 2 
in 0.2mW/mm 2 steps, and the restoring torque for different 
initial twist angle 60; (d) the corresponding paths of power- 
dependent twist angles under different 60 . 




FIG. 3: Schematic of the experimental set-up. The pump and 
probe signals are combined by a 3 dB combiner. The sample 
is positioned in the centre of the waveguide. 1: vector net- 
work analyser; 2: 3dB combiner; 3: rectangular waveguide; 
4: 20 dB attenuator; 5: power amplifier; 6: signal generator; 
7: sample. 



gle, it naturally leads to nonlinear solutions. As shown in 
Fig.[2fd), the power-dependent twist angles under differ- 
ent initial angles demonstrate the evolution from smooth 
nonlinear to bistable response as #0 departs from the an- 
gle of maximum EM torque. In principle, as #0 moves 
further away from the resonance, more noticeable rota- 
tion and hysteresis effects are expected, but higher pump 
power is required (see the case for #0 = 45°). Such evolu- 
tion of the power-dependent nonlinear response can also 
be observed by fixing the initial twist angle but chang- 
ing the pump frequency, as will be demonstrated in the 
experiment below. 

To confirm the feasibility of the proposed nonlinear 
rotatable meta-atoms we carry out a pump-probe mi- 
crowave experiment. To experimentally realise a strong 
nonlinear or even bistable effect, the restoring torque 
from the wire has to be sufficiently small so that the 
structure can be twisted by a large enough angle within 
the maximum available power. We found that rubber is 
a good candidate, since the shear modulus of rubber is 
of the order of 0.2 MPa - 2.4 MPa [24], which is at least 
three orders of magnitude smaller than it is for other 
polymers. 

A schematic of the experimental setup is shown in 
Fig. [3l We have used two separated copper SRRs (inner 
radius r = 3.2 mm, track width 1mm, copper thickness 
0.035 /im and slit width g = 0.2 mm) printed on Rogers 
R4003 substrates (e r = 3.5, loss tangent 0.0027, substrate 
thickness 0.5 mm). The lower SRR is fixed and positioned 
at the centre of a WR229 rectangular waveguide with 
$ = 0°, and the top SRR is suspended with a thin rub- 
ber wire (radius a = 50/im, length d = 20 mm), so that 
it can rotate about the common axis. The two SRRs are 
aligned coaxially, with a face to face distance of 0.75 mm, 
separated by air. The horizontal positions of the SRRs 
are carefully adjusted and the initial twist angle 60 is 
set at around 70°. The mass of the suspended sample 
is 101 mg, which leads to a 6.1% elongation of the wire. 
The Young's modulus is thus estimated as 2.06 MPa, and 
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Waveguide Measurement Array Simulation 
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FIG. 4: Comparisons of resonant frequencies demonstrated 
in a waveguide experiment and numerically calculated for an 
array, (a) and (b) pump at 3.18 GHz; (c) and (d) pump at 
3.21 GHz; (e) and (f) pump at 3.23 GHz. The stable twist 
angles for the array system are shown on the right axes. The 
power in the array is the power incident on each unit cell. 



the shear modulus follows as G ~ 0.69 MPa. The exper- 
imental realization varies slightly from Fig. [TJ however 
the underlying physical mechanism of the nonlinear re- 
sponse due to the dynamic coupling between electromag- 
netic and elastic properties is the same. This equivalence 
is verified in the supplemental material [21]. 

The transmission spectrum measurements are per- 
formed by a vector network analyser (Rohde and Schwarz 
ZVB-20). The CW pump signal is generated by a signal 
generator (HP 8673B) and is further amplified by a power 
amplifier (HP 83020A) before being sent into the waveg- 
uide. Three pump frequencies (3.18, 3.21 and 3.23 GHz) 
are chosen in order to capture the evolution of the non- 
linear response at different distances from the initial res- 
onance. The pump power is increased in 1 dB steps; for 
each step, the sample reaches steady state after 30 sec- 
onds. The mechanical rotation is quite significant and 
can be visually observed in the experiment, thus ruling 
out other possible nonlinear mechanisms such as heating. 

The experimentally observed transmission spectra are 
shown in Fig. S2 of the supplementary file, and the cor- 
responding resonant frequencies are depicted in Fig.H^a), 
(c) and (e), where the predicted evolution from bistability 
to smooth nonlinearity is clearly shown. The initial reso- 
nance (symmetric mode) without pump is located around 
3.256 GHz, and it red-shifts as the pump power increases, 



which indicates that the twist angle is increased. When 
the pump frequency is at the red tail of the resonance, a 
large spectral "jump" (about three times of the resonance 
linewidth) can be observed when the pump power passes 
a certain threshold value [Fig. H(a)]. The thresholds are 
different for increasing and decreasing pump powers. As 
the pump frequency approaches to the initial resonance, 
the spectral "jump" becomes smaller [Fig. HJc)] and fi- 
nally disappears [Fig. Hfe)]. We also observed similar 
effects (not shown) when the pump frequency is at the 
red tail of the antisymmetric mode, in which case the 
two resonances approach each other due to the opposite 
direction of the EM torque. 

To further validate the observed effect, we numerically 
simulated the exact experimental geometry and found 
excellent agreement [2l| . The estimated optimum initial 
twist angle is around 72.5°, and the maximum twist an- 
gles obtained for the three pump frequencies (3.18 GHz, 
3.21 GHz and 3.23 GHz) are around 90°, 85° and 82°, re- 
spectively (U- Finally we remap these angles back to 
the corresponding resonant frequencies according to the 
simulation spectra (see Fig. S3 (b), (d) and (f) of sup- 
plementary file). 

The single meta-atom used in the waveguide exper- 
iment induces image currents in the waveguide walls. 
Thus, our experimental system is analogous to an array 
where neighbouring elements interact. Since the nonlin- 
earity arises from individual rot at able meta- atoms, the 
behaviour in an array will not show qualitative difference 
from a single meta-atom as predicted in the analytical 
model, as long as the neighbour interaction is relatively 
weak. We use full numerical simulation to model an array 
with thickness of a single cell and periodicity of 45 mm 
in the transverse directions, arranged as in Fig. [TJ The 
parameters of the meta-atoms used in the simulation are 
the same as in the experiment. The array and the waveg- 
uide experiment show quite good qualitative agreement 
(see Fig. |D[b), (d) and (f) for the power-dependent reso- 
nance and Fig. SI (b) for the EM torque), thus justifying 
that the nonlinear behaviour observed in the waveguide 
experiment is similar to that of the analogous array sys- 
tem. 

The novel nonlinear "meta-atom" described in this 
work proved to possess a sensitive elastic feedback bring- 
ing nonlinearity to the interaction of EM modes of the 
resonators. The resulting nonlinearity and bistability of 
the response was successfully observed in experiments 
and it turns out that these results can be accurately pre- 
dicted with theoretical modelling. We also note that this 
structure is chiral (except in the high symmetry cases 
of 0° and 180° angle between the rings), it should also 
exhibit nonlinear optical activity, an effect which is rela- 
tively weak in natural media [25| , but can be quite strong 
in metamaterials (26| . 

Although the experimental demonstration in this work 
was performed in the microwave frequency range, the 
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general principle of operation is valid at any frequency 
where a resonant response can be excited in such or sim- 
ilar metamaterials elements, and the way to analyse the 
same phenomena in THz or optical range is conceptually 
the same. 

We believe that this work provides a substantial con- 
tribution to the emerging area of optomechanical and 
magnetoelastic metamaterials, and offers an efficient and 
convenient design for practical applications. 

This work was supported by Australian Research 
Council. The authors are grateful to A. A. Sukhorukov 
for helpful discussions. 
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